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ABSTRACT 

It has been proposed that primordial gas in early dark matter halos, with virial temperatures Tvir > lO'^K, can 
avoid fragmentation and undergo rapid collapse, possibly resulting in a supermassive black hole (SMBH). This 
requires the gas to avoid cooling and to remain at temperatures near T ^ lO'* K. We show that this condition 
can be satisfied in the presence of a sufficiently strong primordial magnetic field, which heats the collapsing gas 
via ambipolar diffusion. If the field has a strength above | B |> 3.6 (comoving) nG, the collapsing gas is kept 
warm (T ^ \(f K) until it reaches the critical density Wcrit ~ lO-'cm"-' at which the roto-vibrational states of 
H2 approach local thermodynamic equilibrium. Ha-cooling then remains inefficient, and the gas temperature 
stays near ^ 10'*K, even as it continues to collapse to higher densities. The critical magnetic field strength 
required to permanently suppress H2-cooling is somewhat higher than upper limit of ^ 2nG from the cosmic 
microwave background (CMB). However, it can be realized in the rare > (2-3)(t regions of the spatially 
fluctuating B-field; these regions contain a sufficient number of halos to account for the z w 6 quasar BHs. 

Subject headings: cosmology:theory - black hole physics - molecular processes - magnetic fields. 



1. INTRODUCTION 

The discovery of very bright quasars, with luminosities 
> 10'*^ erg s"', at redshift z ~ 6 in the Sloan Digital Sky 
Survey (SDSS) suggests that some SMBHs as massive as a 
fewx 10' M0 already existed when the universe was less than 
1 Gyr old (see, e.g.. Fan 2006 for a review). The presence 
of these SMBHs presents a puzzle. Metal-free stars, with 
masses ^100 M0, are expected to form at redshifts as high 
as z > 25 (Abel et al. 2002; Bromm et al. 2002; Yoshida et 
al. 2008), and leave behind remnant BHs with similar masses 
(Heger et al. 2003). However, the natural time-scale, i.e. the 
Eddington time, for growing these seed BHs by > 7 orders 
of magnitude in mass is comparable to the age of the uni- 
verse (e.g. Haiman & Loeb 2001). This makes it difficult 
to reach 10^ M0 without a phase of rapid (super-Eddington) 
accretion, unless a list of optimistic assumptions are made in 
hierarchical merger models, in which multiple seed BHs are 
allowed to grow without interruption, and to combine into a 
single SMBH (Haiman 2004; Yoo & Miralda-Escude 2004; 
Bromley et al. 2004; Shapiro 2005, Volonteri & Rees 2006; 
Li et al. 2007; Tanaka & Haiman 2009). 

An alternative class of explanations involves rapid gas ac- 
cretion or collapse. In this family of models, primordial gas 
collapses rapidly into a SMBH as massive as 10^-10*" Mq 
(Oh & Haiman 2002; Bromm & Loeb 2003; Koushiappas et 
al. 2004; Lodato & Natarajan 2006; Spaans & Silk 2006; 
Begelman et al. 2006; Volonteri et al. 2008; Wise & Abel 
2008; Regan & Haehnelt 2009; Shang et al. 2010), possibly 
by accreting onto a pre-existing smaller seed BH (Volonteri & 
Rees 2005), or going through the intermediate state of a very 
massive star (Bromm & Loeb 2003), a dense stellar cluster 
(Omukai et al. 2008), or a "quasistar" (Begelman et al. 2008). 
These so-called "direct collapse" models involve metal-free 
gas in relatively massive (> 10*^ M0) dark matter halos at red- 
shift z > 10, with virial temperatures T^^^ > 10''K. The gas 



that cools and collapses in these halos must avoid fragmenta- 
tion, shed angular momentum efficiently, and collapse rapidly. 
These conditions are unlikely to be met, unless the gas re- 
mains "warm", i.e. at temperatures Tyir ^ 10"*K. In particular, 
in recent numerical simulations, Shang et al. (2010) found 
that the gas in such halos, when collapsing in isolation, forms 
H2 efficiently, and cools to temperatures of T ^ 300 K. Al- 
though no fragmentation was seen, the gas could ultimately 
fragment on smaller scales that have not yet been resolved 
(e.g. Turk et al. 2009). More importantly, even if fragmenta- 
tion was avoided, the cold gas flows inward at low velocities, 
near the sound speed of ^ 2-3 km s~', with a correspond- 
ingly low accretion rate of ^ 0.01 Mq yr~'. This results in 
conditions nearly identical to those in the cores of lower-mass 
minihalos; extensive ultra-high resolution simulations have 
concluded that the gas then forms a single ^100 M0 star 
(Abel et al. 2002; Bromm et al. 2002; Yoshida et al. 2008) 
or perhaps a massive binary (Turk et al. 2009), rather than a 
supermassive star or BH. 

H2-cooling in early galaxies may be avoided when the gas 
is exposed to an intense UV flux J, either directly photo- 
dissociating H2 (in the Lyman-Werner bands near a pho- 
ton energy of ^ 12eV) or photo-dissociating the intermedi- 
ary H" (at photon energies > 0.76eV). Requiring the photo- 
dissociation timescale, fdiss oc J~^, to be shorter than the H2- 



formation timescale, ffo, 



ex p 



generically yields a critical 



flux that increases linearly with density, y^"' oc p. Since the 
gas in halos with Tyir > lO'^K can cool via atomic Lyman a 
radiation and loose pressure support, it inevitably collapses 
further As a result, in these halos, the critical flux is high, 
y<="' K, \{f-- IQr' (depending on the assumed spectral shape; 
Shang et al. 2010; see also Omukai 2001 and Bromm & Loeb 
2003 who found similar, but somewhat higher values), which 
exceeds the expected level of the cosmic UV background at 
high redshifts. Only a small subset of all T^ir > lO'^K halos. 



which have unusually close and bright neighbors, may see a 
sufficiently high flux (Dijkstra et al. 2008). In order to avoid 
fragmentation, the gas in these halos must also remain essen- 
tially free of any metals and dust (Omukai et al. 2008), which 
may be incompatible with the presence of such nearby lumi- 
nous galaxies. 

In this paper, we consider a different possibility to keep the 
gas warm, relying on heating by a primordial magnetic field. 
Several mechanisms have been proposed to produce a global 
primordial magnetic field (PMF), with a field strength of or- 
der 1 (comoving) nG, during inflation and/or during various 
phase transitions in the early universe (see, e.g. Widrow 2002 
and references therein). If present, the PMF can be strongly 
amplified by flux-freezing inside a collapsing primordial gas, 
affecting H2-formation and cooling. Sethi et al. (2008; here- 
after S08) have shown that 0.2-2 nG fields can significantly 
enhance the H2 fraction during the early stages of collapse. 
Schleicher et al. (2009) found similar results, and empha- 
sized that at the high densities (n > 10^ cm"^) corresponding 
to later stages of the collapse, magnetic heating from 0.1 -InG 
fields results in significantly elevated temperatures. 

In this paper, we consider field strengths higher than those 
in the two previous studies by S08 and Schleicher et al. 
(2009). Our main result is that, in analogy with the UV- 
irradiation, there exists a critical magnetic field, which leads 
to a bifurcation of behaviors. If the PMF has a strength above 
B > Bcrit ~ 3.5 nG, then the collapsing gas is kept warm 
(T ^ lO'* K) until it reaches the critical density n„it ~ lO-'cm"'' 
at which the roto-vibrational states of H2 approach local ther- 
modynamic equilibrium. H2-cooling then remains inefficient, 
and the temperature stays near ^ 10"^K, even as the gas col- 
lapses further On the other hand if B < Bcrit, H2-cooling is 
delayed, but the gas eventually cools down below ^ lOOOK. 
The critical magnetic field strength we find is a factor of 
^two higher than the existing upper limit from the CMB 
anisotropics (Yamazaki 2010; see more discussion in § 5 be- 
low). However, it can be realized in the rare > (2-3)(T re- 
gions of the spatially fluctuating B-field. As we argue, the 
abundance of halos located in these high-field regions is suf- 
ficient to explain the number of the z w 6 quasars observed in 
the SDSS. 

2. CHEMISTRY AND THERMAL EVOLUTION OF COLLAPSING 
PRIMORDIAL GAS 

The thermal and ionization history of gas collapsing into a 
high-redshift halo, in the presence of a primordial magnetic 
field, is described in detail in S08. The dissipation of the mag- 
netic field owing to ambipolar diffusion and to decaying tur- 
bulence in the post recombination era can substantially alter 
the ionization fraction and temperature of the gas, even be- 
ginning before the halo collapse (Sethi & Subramanian 2005 
[hereafter SS05]; Yamazaki et al. 2006; S08; Schleicher et al. 
2009). 

Here we extend the analysis of S08 up to a higher particle 
density in the collapsing halo, and explore higher magnetic 
field strengths. In addition, we made the following changes: 
(a) we track the density evolution of gas in a collapsing halo 
with the model of Dekel and Birnboim (2006; this affects the 
compressional heating rate), (b) we updated the H2 chemistry 
network; specifically, we use the recent compilation in Shang 
et al. (2010; the most significant change is an increase in the 
collisional H2-dissociation rate). 

The evolution of the ionization fraction (x^), magnetic field 
energy density (Eb), temperature (T), and H2 molecule frac- 



tion (XH2) are described by the equations 
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The symbols here have their usual meaning; further details 
and relevant processes are discussed in S08. We list here only 
the processes that have been added or updated. Most of these 
processes relate to the formation and destruction of H2. The 
net rate of formation of H2 through the H" channel is given 
by: 

, kgk\oXHinb 

kioXHini,+k^+(ki2 + k2i)Xpnb+kigXenh + k2oXHinb 

The notation of the reaction rates follows the appendix of 
Shang et al. (2010), except that k^ is the rate of destruction of 
H~ by CMB photons, which, in our case, is important before 
the collapse regime (eq. 8 in S08). The net destruction rate of 
H2, kdes, is: 

^des = kisXHi + kuXp + kiiX^. (6) 

The dominant reaction rates for the range of ionization frac- 
tions, densities, and temperatures we obtain in the entire range 
of our models are: kg, kiQ, km, k\'i, and ^15' 

The cooling processes that dominate Lcooi in primordial gas 
in the density and temperature range we consider are: (a) 
Compton cooling fc,c (eq. 15 in S08), (b) atomic H cooling 
(eq. 16 in S08), and (c) H2 molecular cooling (Galli & Palla 
1998). The heating rate Lheat is given by the magnetic field de- 
cay owing to decaying turbulence {dE^ / dt\urb and ambipolar 
diffusion (liiiB/t/Oambi- In practice, we find that ambipolar 
diffusion always dominates in our case; the dissipation rate is 
given for this process by (Cowhng 1956; Shu 1992; SS05): 
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All quantities in eq. (|7| are expressed at redshift z = 0. The 
time-dependence of the decay rate is given by f{t) = {\+z)'^ 
during the pre-collapse stage, and /(f) c>c p"^/^ during the 
collapse phase. Here M{k) = Ak" is the power spectrum 
of the tangled magnetic field, with a large-^ cut-off at fc = 
A;inax — 235(1 nG/B())(comoving)Mpc"'; fc^^x is determined 
by the effects of damping by radiative viscosity during the 
pre-recombination era (Jedamzik, Katalinic, & Olinto 1998, 
Subramanian & Barrow 1998a; see SS05 for further details). 
Bq, referred to as the magnetic field strength, is defined as 
the r.m.s. value at fc = IMpc"'. Throughout this paper, we 
use n = -2.9 (for justification and further discussion see SS05, 
S08 and references therein). Unless specified otherwise, the 
time-evolution of the magnetic field is assumed to be given by 
flux-freezing, which implies a power-law dependence B oc p" 
on the gas density with a = 2/3. In practice, this scaling may 

' We note here that the rate ki4, assumed to be zero here, and k[g could 
become more important if we used the rates as given by Schleicher et al. 
(2009) and Capitelli et al. (2007). If these rates are used, then the critical 
value of Bq needed to destroy H2 in the collapsing halos decreases to ~ 3 nG 
(cf. Figure 1 below). 



be less steep; below we will explore how our results change 
for different values of a. 

Another important scale is the (comoving) magnetic Jeans 
scale, kj ~ 15(lnG/Bo)Mpc-i (see e.g. Eq. 6 in SS05 and 
references therein). This scale, along with the thermal Jeans 
scale determine the condition that allow a halo to be gravi- 
tationally unstable. In our case, the magnetic Jeans scale is 
larger, and is therefore more restrictive, for Bq > I nG (SOS, 
Figure 4). Primordial magnetic fields can also induce the for- 
mation of first structures in the universe (Wasserman 1978, 
Kim et al. 1996, Gopal and Sethi 2003, Subramanian and 
Barrow 1998a, SS05). The total (dark matter + gas) mass of 
the first structures, determined by the magnetic Jeans length, 
is M ~ 5 X 10*^(50/1 nG)3M0 (e.g. Figure 7 in SS05). 

3. RESULTS AND DISCUSSION 

We show in Figure 1 , 2, and 3 the evolution of the temper- 
ature, the H2 fraction (defined as the ratio of the H2 number 
density to the total hydrogen number density, nHz/wn), and the 
ionized fraction for a single halo from z — 800 (corresponding 
to the initial number density n ~ lOOcm"^ on the left of Fig- 
ures 1, 2, and 3), down to a maximum density of n ~ 10*" cm"-' 
in the collapsed halo. Note that the evolution on these figures 
is monotonically to the right: the x-axis shows the density de- 
creasing to the right (until the turnaround redshift), and then 
increasing again as the halo collapses. 

The figures show the interplay between several physical ef- 
fects. First, the magnetic field decay directly increases the 
temperature. This increases the collisional destruction rate of 
H2, but it also increases the electron fraction (owing to more 
rapid collisional ionization). The larger electron fraction then 
tends to increase the molecular hydrogen fraction, competing 
with the effect of the increased collisional H2-dissociation. 
The molecular hydrogen cooling rate depends on the temper- 
ature directly, and also on the molecular hydrogen fraction. 
As the temperature reaches > 8,000K, atomic cooling domi- 
nates, which, again, is governed by the ionized fraction. 

Figures 1, 2, and 3 should be viewed together to appreci- 
ate the net outcome of these effects on the thermal evolution, 
for different specific values of the magnetic field. A higher 
magnetic field strength generally results in more rapid heat- 
ing, and therefore a higher temperature, at least until halo 
collapse begins. This, in turn results in higher ionized frac- 
tions. However, the thermal evolution is more complicated in 
the collapsing stage, once molecular hydrogen becomes the 
dominant coolant. 

Prior to turnaround, the gas temperature increases roughly 
monotonically with the strength of the magnetic field, up 
to 10'*K. However, since a higher temperature also means a 
higher ionized fraction (Fig. 2), molecular hydrogen forms at 
a faster rate. Once the molecular hydrogen becomes the dom- 
inant coolant, at « ~ 1 cm~^, the halo with a higher magnetic 
field can cool below the less magnetized halo. 

The most interesting case corresponds to Bq = 4nG (shown 
by the thick solid black curves). In this case, the magnetic 
field dissipation rate is high enough to prevent the formation 
of molecular hydrogen in the collapsing halo. As seen in Fig- 
ure 2, the molecular hydrogen starts getting destroyed in the 
collapsing halo in this case. The direct impact of this predica- 
ment on the thermal state of the gas is that the gas fails to cool 
via H2 - the cooling remains dominated by the line excitation 
of the atomic hydrogen. 

Our results are in reasonable agreement with the results of 
Schleicher et al. (2010) for Bq < I nG, the range of magnetic 



field strengths they considered. We note that the molecular 
hydrogen fractions, past the turn-around epoch, are lower by 
more than an order of magnitude compared to the results of 
S08. This is in line with the observation of Schleicher et al. 
(2009) that S08 had underestimated the destruction rate of H2. 
The thermal evolution of the gas in the zero magnetic field 
case (Bo < 10"^ nG in the figures) also agrees with the results 
of the one-zone models studied by Omukai (2001). 

The most important new result of our analysis is the ther- 
mal state of the gas owing to the destruction of the H2 in the 
collapsing halo for B > 4nG. Omukai (2001) analyzed the 
destruction of H2 in the collapsing halo in the presence of 
background UV flux. That analysis suggested that the H2 for- 
mation can be prevented if the halo could be kept at a temper- 
ature ~ 10^ up to a critical density n ~ 2000. At higher den- 
sities, collisions with hydrogen atoms is more effective in de- 
stroying H2 from higher vibration levels, because the relative 
occupation probability of vibration states of the H2 molecular 
approach thermal equilibrium at these densities (the rate A:i5 
in the discussion in section 2, Martin et al. 1996). In the anal- 
ysis of Omukai (2001), the reaction rates are also affected by 
the presence of a background UV field. Our analysis shows 
that the same result (i.e. lack of any H2 cooling) can also be 
obtained by the dissipation of tangled magnetic fields. 

In Figures 4 and 5, we show the heating and cooling rates 
during the collapse regime, for two values of the magnetic 
field, 1 and 4 nG, respectively. Apart from a brief initial pe- 
riod, when atomic cooling roughly balances adiabatic heating, 
the dominant heating and cooling processes in the Bq= I nG 
case, in the collapse stage, are magnetic heating and H2 cool- 
ing, respectively. H2 cooling quickly becomes more important 
than magnetic heating, resulting in a rapid temperature drop 
at densities near n w 0.1-0.2 cm"-' (see Fig. 1). As the halo 
collapses further, the gas begins to recombine and the ion- 
ized fraction decreases (see Fig. 3), the magnetic field dissi- 
pation rate due to ambipolar diffusion increases as oc pjj /p,, 
while the H2 cooling rate grows as Xy,^pl- As a result, the 
magnetic heating catches up with H2 cooling, resulting in a 
nearly constant temperature (see Fig. 1) - however, this oc- 
curs only after the collapse has proceeded beyond the critical 
density n « lO-' cm"-'. Atomic cooling or adiabatic heating 
do not play an important role in the thermal evolution for this 
strength of the magnetic field. 

In contrast, for Bq ~ 4nG, as shown in Figure 5, the mag- 
netic heating roughly balances atomic HI cooling during the 
collapse stage, resulting in a nearly constant temperature T « 
10''K throughout the entire evolution of the contracting gas. 
The magnetic heating gives rise to a higher ionization fraction 
(Fig. 3), and therefore aids the formation of H2. However, the 
high temperature of the gas causes the H2 to get destroyed 
(Fig. 2; this effect is also seen for Bo — 1 nG just before the 
onset of the halo collapse). The net result for this high value 
of the magnetic field is that H2 cannot form fast enough to 
ever become an important coolant. The halo remains at a tem- 
perature ~ lO'' K up to the critical density; as a result, the H2 
fraction is strongly diminished for the subsequent evolution of 
the halo. By experimenting with several intermediate values 
of the magnetic field, we have found that this clear-cut bifur- 
cation in the thermal evolution occurs at a critical magnetic 
field strength of B^it = 3.6 nG. 

As noted in § 2 above, we adopted the model of Dekel and 
Birnboim (2006) to track the density evolution of the collaps- 
ing halo. However, we have checked the robustness of our 



results against a wide range of collapse histories. Specifi- 
cally, after the turn-around stage, we artificially multiplied 
the rate of increase in the density, relative to the Dekel & Bim- 
boim model, by a factor of 0. 1 or 1 . We have found that the 
weaker/stronger adiabatic heating delayed/advanced the onset 
of the catastrophic H2 cooling (e.g. as seen at « « 10~"'^cm"-' 
in Fig.l) to higher/lower densities, but the other qualitative 
features of our results were essentially unchanged. In partic- 
ular, a bifurcation of behaviors was still found, with a critical 
density of Bcrit = 3.6 nG. 

Another uncertainty concerns our assumption of flux- 
freezing. Schleicher et al. (2009) note the possibility of the 
breakdown of this approximation in a collapsing halo. If the 
field is not sufficiently tangled, collapse can occur with little 
dissipation in the direction of the field lines; the magnetic field 
might grow less rapidly than our adopted p^l^ . Our computa- 
tions in the range a ~ 0.55-0.6 show that the critical magnetic 
field required to prevent the halo from cooling increases to 
Bo — 5-7 nG. This also leads to the interesting possibility that 
the average Bo — a fewnG, required to form SMBHs, might 
be detectable by the future CMB experiments (e.g. Yamazaki 
et al. 2010) and the 4-5(t fluctuations of the field might leave 
their trace in the formation of SMBH. We hope to explore this 
possibility in the future. 

At present, the best upper limits on the PMF are in the range 
2-3 nG from CMB temperature and polarization anisotropics 
and from early structure formation (Subramanian & Barrow 
1998b, 2002, Durrer, Ferreira & Kahniashvili 2000, Seshadri 
& Subramanian 2001, Mack, Kahniashvili & Kosowsky 2002, 
Lewis 2004, Gopal & Sethi 2005, Kahniashvili & Ratra 2005, 
Giovannini & Kunze 2008, Yamazaki et al. 2008, Finelli, Paci 
& Paoletti 2008). Recently, Yamazaki et al. (2010) obtained 
bounds on the strength and spectral index of the magnetic field 
power spectrum: Bq < 2nG and n < -1 .4. Our results cannot 
be compared directly to their analysis as we use only a sin- 
gle value of spectral index (n = -2.9) here. The upper limits 
from CMB on Bo for this value of n are considerably weaker 
(Lewis 2004). We note that the effect we discuss here might 
be pronounced for a larger n for a given Bo. More detailed 
analysis would be required to directly compare our results 
with Yamazaki et al. (2010). Our analysis suggests that for 
B() > 4nG, the fragmentation of collapsing halos could be pre- 
vented, and therefore this could be considered an independent 
upper bound on the value of Bq- 

4. THE MASS OF THE CENTRAL OBJECT 

The mass of the central object that forms in a collapsing 
proto-galaxy can be approximated as follows: there exists 
a radius at which the mass accretion time-scale face equals 
the Kelvin-Helmholtz time scale Tkh for a proto-star, with 
the proto-stellar mass equal to the gas mass enclosed within 
this radius. For metal-free gas, fxH is approximately 10^ 
years, with only a mild dependence on the proto-stellar mass 
(Schaerer 2002). The expected mass of the central object 
then scales approximately as M ex t~^^ ex c, ex T^^^ (see, e.g., 
Shang et al. 2010 for the last scalings with the sound speed 
and gas temperature). This implies that a stellar mass of ^ 
2OOM0, expected for T = 300K, can increase to « 4 x 10"* M© 
when H2-cooling is inefficient and T « lO'^K (in their three- 
dimensional simulations, Shang et al. find a somewhat still 
steeper scaling). Our proposal here is that a small fraction 
of halos at z = 10- 15, which contain pristine, metal-free gas 
when they collapse, and which reside in regions of an unusu- 
ally high initial seed magnetic field, may produce a SMBH 



with a mass of up to ^ 10"*"^ M0. The time available between 
z = 6 and z = 10- 15 is w (4-6) x 10^ yrs, allowing for a fur- 
ther growth in mass by a factor of w (2 x 10^) -(3 x 10*") at the 
e-folding time of 4 x 10^ yr, (corresponding to Eddington- 
limited growth at the radiative efficiency 10%). Hence, the 
10''"^ M0 BHs, produced through the primordial magnetic 
field, can indeed grow into the > 10^ M© SMBHs by z = 6. 

The smallest total (dark matter-ngas) mass that can col- 
lapse for Bo ~ 4nG is M ~ 3 x lO'^M© (as mentioned in 
section 3 above). These halos could form either as a result 
of gravitational instability in the standard ACDM model, or 
via PMF-induced density perturbations (for details on the lat- 
ter scenario, see e.g. SB08). The abundance of halos in the 
PMF-induced structure formation case drops very sharply for 
masses above the Jeans mass (e.g. SS05), and for simplic- 
ity, we conservatively drop this contribution in our analysis. 
In the usual ACDM model, using the fitting formula for the 
halo mass function from Jenkins et al. (2001), and the current 
best-fit cosmological parameters from Komatsu et al. (2009), 
we find that the abundance of all M > 3 x 10'° M© halos at 
z=10is«5x 10"^ (comoving) Mpc~ . At the somewhat 
higher redshift of z = 15, the abundance of halos above the 
same mass drops sharply to « 3 x 10"^ cMpc"^. 

The space density of > 10^ M© SMBHs, inferred from 

the observed abundance of bright z ~ 6 quasars, is ~ Cq 
(comoving) Gpc"^. Here eg denotes the duty cycle, defined 
as the fraction of the Hubble time that z = 6 supermassive 
black holes are observable as luminous quasars. Assuming 
a quasar lifetime of ^ 50 Myr (e.g. Martini 2004), we have 
eg ^ 0.05, and the space density of z = 6 SMBHs is ^ 20 
Gpc~^. Therefore, at redshift z = 10, a fraction as low as 
/ - (20 Gpc"^)/(5 X 10"^ Mpc"^) = 3 X 10"^ of the whole pop- 
ulation of M > 3 X 10'" M0 halos is sufficient to account for 
the presence of these rare z ~ 6 SMBHs. Such a small frac- 
tion could would correspond to ^ 2.8cr upward fluctuations of 
a Gaussian random PMF. Because of the strong reduction of 
the relatively high-mass halos at higher redshift, essentially 
every M > 3 x 10^° M© halo would have to host SMBHs at 
z « 15, to match the comoving abundance of z = 6 quasar BHs. 
z ~ 15 is therefore the earliest epoch for forming the heavy 
SMBH seeds as envisioned here. As long as the PMF ampli- 
tude is B > 1 .2nG, the critical flux we find, B « 3.6nG, would 
be reached by the ^ 3cr upward fluctuations. 

As mentioned above, our results are analogous to earlier 
work, proposing that a UV flux keeps the gas hot, and results 
in SMBH formation. However, it is worth emphasizing two 
important differences in these two scenarios. First, as Figure 1 
shows, even though there is a clear critical B-field value, di- 
viding the thermal evolution into two different regimes, even 
for B-field strengths below this value, the minimum tempera- 
ture reached by the gas is significantly elevated. Hence, sub- 
critical magnetic fields can still significantly increase the mass 
of the central object. Second, the gas in the halos must remain 
essentially free of any metals and dust, in order to avoid frag- 
mentation. This may well be impossible to reconcile with the 
large required flux, which has to arise from a close neighbor 
galaxy, in the UV-irradiation scenario. On the other hand, 
in the magnetic-heating case, there is a-priori reason why a 
region where the PMF has an upward fluctuation would be 
more likely to be metal enriched. In fact, quite the contrary: 
the large magnetic Jeans mass will suppress the collapse of 
gas into lower-mass halos at high redshift, naturally exclud- 



ing any prior metal-pollution farther up in the merger tree of 
the halo in which the putative SMBH forms. 

5. CONCLUSIONS 

We found a plausible novel mechanism to form high- 
redshift SMBHs by direct gas collapse in early dark matter 
halos, aided by heating from the dissipation of a primordial 
magnetic field. The model avoids many of the assumptions 
required in earlier models (such as an extremely high UV flux 
and the absence of H2 and of other molecules and metals), 
but it does require a large primordial magnetic field, and re- 
lies on metal-free primordial gas. We expect that, in gen- 
eral, any other heating mechanism, which can compete with 
atomic HI cooling in the collapsing halo, down to a density of 
n ^ 10^ cm"^, would produce the same effect as the B-field 
utilized here. 

Interestingly, our model requires a magnetic field strength 
that is close to the existing observational upper limits. We 



therefore expect that the values Bo — a few x 10"^ G, required 
to form SMBHs, could be detectable by the future CMB ex- 
periments. The upward fluctuations of such a strong primor- 
dial magnetic field might also leave their trace on cosmic 
structure formation: because of the high magnetic Jeans mass 
in these regions, the formation of dwarf-galaxy sized halos, 
with masses below M ~ a few x IQ^^^Mq, will be prevented. 
We expect this can lead to further constraints on our SMBH- 
formation scenario, analogous to constraints on warm dark 
matter models (which produce a similar suppression of small- 
scale structures; Barkana, Haiman & Ostriker 2001). 
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Fig. 1 . — The temperature evolution of a patch of the intergalactic medium is shown as it initially expands and then turns around and collapses to high density. 
The different curves correspond to different values of the assumed primordial magnetic field, as labeled. The gas evolves from the left to the right on this figure. 
The left panel shows the expanding phase, starting from an initial density of Si 100 cm"' (corresponding to the mean density at redshift z — 800) and ending at 
the turnaround just below n = 10"^ cm"'. The right panel follows the subsequent temperature evolution in the collapsing phase. 
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Fig. 2. — The evolution of the H2 fraction in the same gas clouds shown in FigurefT] 
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Fig. 3. — The evolution of the ionized fraction in the same gas clouds shown in FigurefT] 
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Fig. 4. — The heating and cooling rates are shown for various processes as labeled, forB = 1 nG. The rates are in the units dt/d7, = Hq'Q„ (l+z) ^"- 
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Fig. 5. — Same as Figure 4, but for B = 4 nG. 



